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Abstract. We consider the two-point function of nucleon current in nuclear matter and write a QCD sum 
rule to analyse the residue of the nucleon pole as a function of nuclear density. The nucleon self-energy 
needed for the sum rule is taken as input from calculations using phenomenological AW potential. Our 
result shows a decrease in the residue with increasing nuclear density, as is known to be the case with 
similar quantities. 
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1 Introduction 



average: For any operator O, 



An important topic in strong interaction at non-zero tem- 
perature and chemical potential is the propagation of hadrons 



Z ^ Tr 



-l}(H-tiN) 



through these media. Their modified couplings and self- 
energies are useful not only in analysing the experimental 
data on heavy-ion collisions, but also in extracting indica- 
tions of an eventual phase transition in the medium. 

At low density, such a propagation in the hadronic 
phase can, in principle, be studied by invoking chiral per- 
turbation theory to evaluate the appropriate two-point 
function [HI]. This effective theory of QCD is of par- 
ticular advantage at low temperature (and zero chemi- 
cal potential), when the heat bath is dominated by pions. 
Their couplings with themselves and with other hadrons 
are highly constrained by the chiral symmetry of (mass- 
less) QCD. Thus not only are the vertices with pions re- 
lated to those without the pions, but also they are sup- 
pressed by powers of pion momenta. Using these vertices, 
one can get the hadron parameters to a good accuracy 
by evaluating only a few relevant Feynman graphs. It has 
been applied with much success in calculating the proper- 
ties of the pion and the nucleon at finite temperature [31 

SIEHZIIH]. 

In this work we are interested in finding the effect of 
strong interaction on the propagation of nucleon at non- 
zero nucleon chemical potential (/x) i.e. in nuclear matter, 
at zero temperature {(3~^ = 0). Quite generally, this prop- 
agation is studied by considering the two-point function, 

n{E,p) =i I dtSxe'^'^*-P-''\Tri{x)f]{Q)) (1) 



where H and N are the Hamiltonian and the number op- 
erator of the system. 

Unfortunately, a straightforward calculation of Feyn- 
man graphs in this case, similar to the one at finite tem- 
perature, is not possible at present. The difficulty is due 
to the appearance of new and presumably large couplings. 
They are shown in Fig. 1, where the complete 77-nucleon 
vertex and the complete nucleon self-energy are analysed 
in terms of low order perturbative vertices, whose struc- 
tures may be obtained by the methods of chiral perturba- 
tion theory. 

Thus, up to terms proportional to a single nucleon field 
il){x), the current ri{x) is given by |11| 



77(x) = A 1 + 



i4>{x) 



2F^ 



-75 



ipix) 



(2) 



where </>(a;) is the pion field and the dots stand for terms 
with increasing number of pion fields. Here F^^, the so- 
called pion decay constant, is defined by the vacuum-to- 
pion matrix element of the axial- vector current A^^{x), 



F^=93MeV, (3) 



just as A is defined by the vacuum-to- nucleon matrix ele- 
ment of the nucleon current ri{x) 



{0\ri{x)\N{p)) = Xu(j))e 



ip-x 



(4) 



of a three-quark current ri(x), having the quantum num- 
bers of the nucleon pillOj. Here (• • •) denotes ensemble 



where u{p) is a positive energy Dirac spinor. The value of 
A is obtained from QCD sum rules for nucleon in vacuum 
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(1.2 ±0.6) X lO-^GeV^ 



(5) 
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Fig. 2. Graphs representing the two-point function in the low 

energy region. 



Fig. 1. Complete vertex and complete self-energy (denoted by 
solid circles) analysed in terms of loop graphs. 



Terms in 77 proportional to 'ijjijjtp may also be obtained in 
the same way, bringing in two more new coupling con- 
stants [TT]. We now see in Fig. 1 that unlike the vertex 
rjipcj) in graph (b), which is related to the vertex rjip of 
graph (a) itself, the vertex rjipipil) in graph (c) is unknown 
and unlikely to be small. 

For the self-energy graphs, we note the pion-nuclcon 
interaction Lagrangian [12j . 



C^nt = -^-0(x)7p75rXa;)a^(/)*(a;) 



(6) 



where qa is the axial-vector coupling constant of the nu- 
cleon, gA = 1.26. We again see in Fig. 1 that unlike 
the self-energy graph (e), which can be calculated with 
this Lagrangian, the other graph (f) poses difficulty, even 
though chiral symmetry dictates the form of the four- 
nuclcon effective Lagrangian, whose coefficients can be 
fixed by the AW scattering lengths jl3] . Indeed, if one does 
calculate the graph (f) with this four-nucleon interaction, 
one gets an unacceptably large value for the nucleon self- 
energy [Mj. The problem can be traced to the fact that 
there are bound and virtual states very close to threshold 
in the NN system. 

In view of these difficulties, we give up calculating the 
nucleon self-energy and content ourselves with evaluating 
only the nucleon pole residue, by writing a QCD sum 
rule for an appropriate amplitude representing the two- 
point function. (There are several approaches to QCD 
sum rules in medium [l5lfT6lfT7] . The one closest to ours 
is of Ref.[T7].) To determine the nucleon pole term in the 
medium, we still need the nucleon self-energy, which we 
take from the variational and the Brueckner type calcula- 
tions using the phenomenological NN interaction poten- 
tials [181119] . The contributions from the remaining low 
energy singularities (branch cuts) are obtained by evalu- 
ating the relevant Feynman graphs. 

Such sum rules in medium arc usually written for com- 
plete amplitudes, including their contributions from the 
vacuum. As noted in Ref. [20j . a sum rule of this type can- 
not, however, determine the pole residue quantitatively. 
The reason is that in such sum rules the vacuum contri- 
butions dominate over the density dependent ones, at least 
at low densities, making the latter appear as non-leading 
terms. 

In this work we subtract out the vacuum sum rule 
from the corresponding one in medium, that is, we ex- 
clude the vacuum contributions from both the spectral 
and the operator sides. While this exclusion of the leading 
terms is expected to make the sum rule more sensitive to 



the medium dependent quantities, it also calls for a more 
careful scheme to saturate it. 

We construct the spectral and the operator sides of the 
sum rule in Sec. 2 and 3 respectively and write the sum 
rule for the nucleon pole residue in Sec. 4. It is evaluated 
in Sec. 5. Our discussions arc contained in Sec. 6. 



2 Spectral side of sum rule 

The low energy part of the spectral side of the sum rule is 
constructed from singularities generated by the low mass 
physical particles. Besides the complete nucleon pole, we 
also include branch cuts from ttN exchange. We are thus 
led to consider the contributions of graphs of Fig. 2. As 
already stated, the solid circles in graphs (a), (b) and (c) 
are not for diagrammatic evaluation: The one represent- 
ing the current-nucleon coupling is left as unknown, to be 
determined by the sum rule. The other standing for the 
self-energy is taken from the results in Ref. |18[[T9] . The 
remaining graphs do not involve any new couplings. Con- 
cerning graphs (e) and (f) we take only their non-pole 
parts, as the pole part is already included in graphs (a) 
and (b). 

We shall calculate amplitudes in the real time version 
of the field theory in medium, where a two-point function 
assumes the form of a 2 x 2 matrix [21] . (This formulation 
is reviewed in |22j). But the dynamics is given essentially 
by a single analytic function, obtained by diagonalising the 
original matrix. Thus if 7Iii(£',p) is the 11-component of 
the original matrix amplitude, the corresponding analytic 
function, to be denoted by the same symbol as in Eq. (1), 
has the spectral representation [23], 



n{E,p)^ ( 

J — c 



dE' a{E',p) 
2tt E' -E- i7je{E') 



(7) 



where the spectral function is related to the imaginary 
part of Till by 

a{E,p) = 2coth{P{E - fi)/2}lmnii{E,p). (8) 

For generality we calculate the amplitudes retaining both 
fi and P and take the limit of zero temperature later. 

The two-point function due to the free propagation of 
nucleon, namely 

A2 
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is modified by the vertex and the self-energy corrections 
of Figs. 2 (a), (b) and (c) to 



niE,p)\,, 



X* 



{a+b+c) 



(9) 



where A* is the modified couphng parameter in nuclear 
matter and S{p) is the nucleon self-energy acquired in 
this medium. 

In this work we restrict to p = 0, when S has the sim- 
ple Dirac decomposition, S — Us + 7°^v- Then Eq. (9) 
may be rewritten as 



n{E)\i^a+b+c) - -A* 



{E - m{){E - 1112) 



(10) 



The scalar part Es of the self-energy changes the mass m 
of the free particle to the effective mass m* in the medium, 
m* = m + Us- Similarly the vector part Sy shifts the rest 
energies, ±m of the free nucleon and the antinucleon, to 
m-i = m* -t- Sv and TO2 = — m* -|- respectively of the 
corresponding quasi-particles. Following our discussion in 
the Introduction, we work with the subtracted nucleon pole 
term, 



n{E)\^a+b+c)=n{E)\(, 



+b+c) 



+ A 



E^ 



(11) 



Here and below we use a bar over amplitudes U and spec- 
tral densities a to denote subtraction of the corresponding 
vacuum contributions. 

We now evaluate the remaining graphs of Fig. (2). The 
imaginary part of the 11 component of the matrix ampli- 
tude for graph (d) is given by [23,24,25 

Imn(£:)ii|(rf) ^ /SA^ttA r d^q 
tanh[/3(£; - /i)/2] ~ 1^ / J (27r)34coiW2 
[(— 7°tJi + m){{l - n+ + n)S{E - UJ\- uj-i) 
+ (n+ + n)&{E - tJi + 0)2)} - {^1.2 -^1.2, n+ 

(12) 



where wi = i/m^ + q^^ ^2 = ^y rri^ + q'^ and n± and 
n are respectively the distribution functions for nucleons, 
antinucleons and pions, 



n±{u;i) 



1 



n{LU2) = 



1 



(13) 



Terms without the n's are the vacuum contributions, which 
we subtract out, as already stated above. Further, we re- 
strict to zero temperature so that we have to calculate 
only the term proportional to n+ — > 0(/i — lui): we thus 
get 



3A2 



16ttF2E 



(14) 



on the Landau cut and the negative of the same quantity 
on the unitary cut. Here uj is the nucleon energy expressed 
in terms of the total energy E, uj = {E"^ + m? — rn^) / {2E) . 



Because of the ^-function in the integrand, the two cuts 
originally over the regions Q < E < m — rriT^ (Landau) and 
m -\- rriTr < E < 00 (unitary) shrink respectively to 



H — \J p?- — w?- -I- m\ < E < m 



(15) 



and 



m + rriTT < E < fi + -^//i^ — + . (16) 

Then the desired analytic function in the form of Eq. (7) 
is given by 

dE' a(S')|(d) 



n{E)\ 



id) 



C 2n E' -E- irie{E') 



(17) 



where the subscript C denotes the difference of two inte- 
grals, 

(18) 



Landa 



unitary 



Similarly the graphs (c) and (f) give 



n{E)\ 



(e+f) 



2{-f°E - 



£2 - J7l2 



dE' ^{E')\u 



where 



c 2tt E' -E- i7j€{E') 

(19) 



aiE)\ 



^Tj^L ^""' - mmm^-Eu;)+^°m{E-L,)} 

(20) 

Of course, we have yet to subtract out from Eq. (19) its 
nucleon pole part, as the complete nucleon pole contribu- 
tion is already represented by Eq. pT|) . 

At this point we adopt an improved choice of the am- 
plitude following the suggestion in Ref. [5D]. Of the two 
quasi-particle poles, the one corresponding to the antin- 
ucleon cannot be treated as having narrow width, as we 
do in Eqs. (10) and (11)), because it can be annihilated 
strongly in nuclear medium. This ill-represented contri- 
bution may be suppressed by working with a modified 
amplitude. We split U into even and odd parts. 



n{E) 



n'^'\E^) 



and deal with the combination, 



^n^'\E^) 



■ En^''\E'), 



m2n^"\E^ 



(21) 



(22) 



(Recall that m2 is the pole position of the quasi anti- 
nucleon.) Indeed, if we separate the unsubtracted pole am- 
plitude n{E) given by Eq. (10) into even and odd parts 
and combine them as in Eq. (22), we get the amplitude 



2A* 



1 



E^ -ml2 



(1 + 7°) 



where the pole at E' = m2 is removed to £^ = mi. Note also 
that it is proportional to i(l -I- 7"). As we arc interested 
in a sum rule for A* , we shall project all amplitudes on to 
this combination in Dirac space 



^ The amplitudes proportional to |(1 — 7°) are expected to 
be relatively small and so comparable to those left out in our 
approximation. Thus the corresponding sum rule may not be 
reliable, besides the fact that it will not involve A*. 
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Wc thus get the amphtudes corresponding to the dif- 
ferent graphs of Fig. (2) as 

)\{a+b+c) - --f^, — — + —f^^ — zzr- y^-^) 



\{d) 



3A2 



327r2F2 



TT JC 



ml E^- 
dE' f{E') 
E' E'^ - E^ 



(24) 



3AV f dE'E' + m f{E') 



16tt^F^ Jc E' E' E'^ - E'^ 



Here f{E) is given by 



f{E) = {E — m2)\/ LO^ — ni^{uj — m) 



(25) 



(26) 



In the last amplitude we have removed the nucleon pole 
contribution from graphs (e) and (f). The sum of these 
amplitudes will give the low energy part of the spectral 
side of the sum rule. 



3 Operator side of sum rule 

In obtaining the operator product expansion, we need the 
explicit form of the quark current ri{x) da with spin and 
isospin indices D and i. Of the two independent possibil- 
ities involving three quark fields, we choose here the one, 
which for proton {i ~ 1) is0, 

where C is the charge conjugation matrix and a, b, c are 
the colour indices. Because of subtraction of the vacuum 
amplitude, the unit operator in the operator product does 
not appear in n{E). As is well-known, at higher dimen- 
sions, there are two sets of contributing operators in the 
in- medium sum rule: the old set, appearing already to the 
vacuum sum sum rule and the new set, involving w^, the 
four-velocity of the medium. From the ensemble average 



^ Let us discuss here the preference, if any, for this choice. 
The current written above corresponds to an axial-vector di- 
quark coupled to a quark. Similarly, the other possibility, where 
7^'s are replaced by a^^'s, corresponds to an axial-tensor di- 
quark coupled to a quark. Reliability of results from QCD sum 
rules requires that the nucleon pole term dominate the contin- 
uum contribution and at the same time the higher order power 
corrections be small compared to the leading ones retained in 
the sum rule. On the basis of these criteria, one may choose the 
current with 7p's [26]. There is a different line of investigation 
of the three-quark currents based on models with instanton 
dominated (non-perturbative) vacuum. For a review see Ref. 
|27) . Each of the two currents mentioned above may be writ- 
ten as a linear combination of two other three-quark currents, 
with relative coejficients ±1, where they couple a scalar and a 
pseudo-scalar diquark with a quark. These models find a large 
coupling of the nucleon with the current containing the scalar 
diquark compared to the one containing the pseudoscalar di- 
quark. Thus they tend to support the use of any of the original 
forms (with and (Tfj,y) in writing the QCD sum rules. 



of the old set of operators, we have to subtract out the 
coirtribution of the vacuum state. 

We denote by q any of the u and d quark flavours. Then 
the contributing operators of lowest dimension are qq, qt^q 
{=q''q in the rest frame of matter, where u^^ = (1,0)). 
Next, there are the operators of dimension four, namely 
G2 = {as/TT)Gl^G^""' and O^^s = ut^u^Ol'S, where the 
(traceless) energy-momentum tensors of quarks and glu- 
ons are given respectively by 



^1^^= i^lli^^'^l- ^9t..^qqj , (27) 

q—u,d 

1 ^28) 



Here m is the average quark mass of u and d quarks and 
GJJj, arc the gluon field strengths. Of the remaining oper- 
ators we retain only the four-quark operators. Below we 
shall work in the limit of SU(2) fiavour symmetry; thus 
(OluitlO) = {Q\dd\0) , (uu) = {dd) , etc. 

In terms of the above operators, the operator expan- 
sion gives 

n{E, 0) ^((^^) + Hu''u)j")E^ H-Eyfi'^) 

-6^(^(G^) + 5(^?0)/£;in(-£;VM^) 



-|^(7°(S«)^ + 2(S«)(^^t^,)) 



(29) 



with coefficients to zeroth order in 0. The renormaliza- 
tion scale n is taken at 1 GeV. 

The four-quark condensates in medium encountered 
above have been factorised in a manner similar to that 
by saturating the vacuum condensates with the vacuum 
intermediate state - an approximation that may not be 
as good as in vacuum. To rectify the error, we adopt the 
suggestion of Ref. [50] to interpolate the factorised (uu)^ 
between its values in vacuum and in medium. Thus we 
shall replace it as 



(1- f){0\uu\Of + f{uuy 



(30) 



where / is a real parameter in the range < / < 1. 

Recalling that the combination (22) of amplitudes is a 
function of E^ , we can go to the spacelike region by set- 
ting E^ = -Q^, Q2 > 0. While ^J.^ = 1 GeV^ is a natural 
scale for the expectation values of the operators, it is not 
convenient for the coefficients, as they contain powers of 
\n{Q^/pi^) in higher orders. As is well-known in the con- 
text of deep inelastic scattering, it is possible to use the 
renormalization group equation to get rid of these (large) 
logarithms by shifting the scale for the coefficients from 

to Q^. The process brings in the (small) anomalous 
dimensions of the operators; also the operators 0"^ and 
0^ will mix PSII^ . However, as we shall vary the Borel 
mass in the neighbourhood of 1 GeV, the renormalization 



^ The lowest order coefficient of is proportional to 
asln(-SV/i') 
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group improvement will give small corrections, which we 
shall ignore in the present analysis. 

The nucleon number density n is related to the Fermi 
momentum pp hyn = 2p|,/(37r^). In normal nuclear mat- 
ter, it is given by no = (110 MeV)'^ corresponding to 
Pf ~ 270 MeV. To first order in n, the change in the ex- 
pectation value of an operator O in nuclear matter relative 
to that in vacuum is given by its nucleon matrix element 
as 

(O) = (0|O|0) + ^P^n (31) 

We apply this equation to the different operators. For uu 
and u^u, we get 

{uu) = (0|uw|0) + ^n, {u''u) = (32) 



where a is the so-called nucleon cr-tcrm, |30| . 

cr = rh{p\uu\p) / m = 45 ± 8McV. 



(33) 



The quark mass and the vacuum condensate are related 
by the Gell-Mann, Oakes and Renncr formula [5T] . 



F^ml = -2m(0|uM|0) . 



(34) 



Two determinations of these quantities exist in the liter- 
ature, namely 



rh = 7.2 McV, 
m = 5.5 MeV, 



|ww|0) 
|mm|0) 



-(225McV)3 (35) 
-(245MeV)3 (36) 



obtained respectively in Refs. [31] and Let us note 
here that using the formula (34), we may write the con- 
densate (32) in nuclear matter as 



(37) 



{uu) = (0|i-m|0) 1 



which vanishes at n = 2.8no. 

Next, for 0^ we can write the nucleon matrix element 



as 



{p\0iM = 2Af{p^p, ~ g^,m^/A) , (38) 

where the coefficient is determined by the first mo- 
ment sum rule for the quark distribution function in deep 
inelastic scattering [34] . Evaluated at the momentum scale 
of 1 GeV, it has the value Af = 0.62 Then noting the 
normalization condition, (0|(9^j^|0) = 0, we get from Eqs. 
(31) and (38), 



{e^) = -mAfn 
4 



(39) 



Finally for the operator G^, wc use the trace anomaly 
to relate it to the trace of the full energy momentum tensor 

O't,^ -l^C^ + 2rhuu + c-l (40) 
^8 

where we add the c-number term to fix again its vacuum 
normalization, (O|0^i/|O) = 0. Taking the vacuum and the 
ensemble expectation values, we get 



(G2) = (0|G2|0)--(m-a)n 



(41) 



With the above results, wc can subtract out the vac- 
uum contributions from Eq.(29) for n{E, 0) and write the 
result for the amplitude combination (22) as 



7T(Q 



2n ope 



A 
"8^' 



37712 , 

^ in 



2Gm2 
~3Q^ 



(42) 



where A, B and G stand for the constants. 



A: 

B 
C 



a 

771 



12 



bmA^ — g(™ ^ f ) 
{Q\uu\Q) f/-^ + 3 

V 771 



(43) 



4 Sum rule 

It is now simple to take the Borcl transform of the spectral 
and the operator sides and get the desired sum rule 



c 



"1 - TO2 ^-,nVA/^ 

2to* 



dE ^, , { f E + TTi 

TO2 



E 



647r2F2m* 



2A2r 



M2 
8^ 



AV. 



87r2 



2Gm2 
3A/2 



(44) 



where f{E) is given by Eq. (26) and 

^1 = 1- e-'^'/*^', ^2 = 1 - (1 + iyVAf2)e-^VM^ 

The deviation of Vi^2 from unity represents the contri- 
bution from the high energy region on the spectral side, 
obtained by continuing the result for operator expansion 
to the time-like region. Here is a parameter determin- 
ing the onset of this continuum contribution. We take 
W — 2 GeV, as assumed for the vacuum sum rules [S]. 

Let us recall here that a special feature of the present 
QCD sum rule is its sensitivity to the medium dependent 
quantities, requiring a more careful saturation of its spec- 
tral and operator sides. Since wc consider nuclear matter 
at zero temperature, it consists only of nucleons. Then, 
on the spectral side, the two-particle density dependent 
contributions can arise from intermediate states consisting 
only of N with p,uj, (j), ■ ■ ■ . As seen from the inequalities 
(15,16), the range of energy over which these states may 
contribute is highly restricted. Thus at normal nuclear 
density, these ranges are only 1080 MeV < E < 1290 MeV 
on the unitary cut and 670 MeV < < 800 MeV on the 
Landau cut. Inclusion of more particles in the intermedi- 
ate state will, of course, increase these ranges, but their 
contributions will fall off exponentially in the integrand. 
So we include only the continuum contribution from the 
ttN intermediate state, in addition to the modified nucleon 
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0.8 1 1.2 1.4 1. 

M 



Fig. 3. Borel mass dependence of A*^/A^ for four relative den- 
sities (n/no), namely 0.25, 0.5, 1.0 and 1.5. The curves spotted 
without and with circles correspond to — .0012 GeV® and 
A^ = .0018 GeV*^ respectively. 

pole. Although the continuum contributions over the Lan- 
dau and the unitary cuts are individually rather big, the 
two largely cancel out in the difference (sec Eq. (18)), leav- 
ing the nucleon pole contribution to dominate the spectral 
side. 

The operator side is also well saturated with the usual 
low dimension operators. The leading operators, qq and 
q^q bring in a large contribution, as represented by the 
term with A in Eqs. (42, 44), compared to those of dimen- 
sions four and six retained in the sum rule. Also the nu- 
merous operators of dimension five and higher do not ap- 
pear to contribute significantly, at least individually 1201 . 

5 Evaluation 

In presenting numerical results, we observe that among all 
the parameters, it is the A, which enters most sensitively in 
the sum rule and also suffers from the largest uncertainty 
in its value (Eq. (5)). Under these circumstances, we fix A 
by requiring maximal stability of the results against vari- 
ation with respect to the Borel mass in a range similar 
to that required for the vacuum sum rules [9]. Numeri- 
cal evaluation shows that bigger the value of A, the more 
stable is the result. We show this stability in Fig. 3 for 
different relative densities, taking A^ = .0012 GeV^, the 
central value and ~ .0018 GeV®, the largest in the al- 
lowed range. (Here we take the self-energies from Ref. [T^; 
those from Ref. |TH] give even better stability of the re- 
sults.) With the latter value of A, we see that there is a 
reasonable plateau up to about normal nuclear density, 
vindicating our saturation scheme. At higher densities, 
however, the plateau seems to disappear, indicating in- 
sufficient saturation. 

We next consider the uncertainties in the results from 
the remaining inputs. First consider the self-energies. As 



Fig. 4. Plot of \*^ /}? as a function of n/no for three Borel 
masses (M), namely 0.8, 1.0 and 1.2 GeV. The upper and 
the lower panel show the results using self-energies given in 
Refs. [18) and [19] respectively. The self-energies are available 
at discrete values of densities, the minimum of which marks 
the beginning of our curves here and in Fig. 5 below. 

already stated, these have been determined by two groups 
[Hmj, using entirely different methods based on phe- 
nomenological NN potentials. Fig. 4 depicts results using 
these two sets of self-energies. Clearly it is not a source 
of any significant uncertainty in the results, at least up to 
normal nuclear density. 

Finally we vary the values for the sigma term and for 
the pair, the quark mass and the quark condensate. As 
seen from Fig. 5, the uncertainty in these parameters again 
does not give rise to any significant spread in the values of 
A*. Also the term with the parameter / in Eq. (43) arising 
from the approximation to the four-quark condensate is 
relatively too small to change the results appreciably. We 
thus show unambiguously a decreasing trend for A* with 
the rise of density at least up to normal nuclear density. 



6 Discussion 

We describe here a method to find A*, the parameter cou- 
pling the nucleon current to the nucleon state in nuclear 
matter. The method consists of writing down a QCD sum 
rule in this medium for an appropriate combination of 
amplitudes representing the ensemble averaged two-point 
function of nucleon current. With the present scheme of 
saturating the sum rule, the results are reliable up to nor- 
mal nuclear density, within which one finds definitely a 
decrease of A* with increase in density. For quantitative 
results beyond this density, we have to improve upon our 
saturation scheme. A test of this improvement at a higher 
density is offered by observing a plateau in A* at such 
density as a function of the Borel mass. If, however, we 
continue the present calculation beyond its range of valid- 
ity, it goes to zero at about twice this density. 



S. Mallik, Sourav Sarkar: QCD sum rule for nucleon in nuclear matter 



7 



1.0 



0.8 



0.6 



5< 



0.4 



0.2 



set I 
set II 
set III 
set IV 



0.0 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

ii/n„ 

Fig. 5. Plot of A*^/A^ as a function of n/no with four different 
parameter sets for a and the pair, rh and {0|?iu|0). I: cr = 
45MeV, Eq.(35); II: cr = 45MeV, Eq.(36): III: a = 53MeV, 
Eq.(35): IV: a = 37MeV, Eq.(35). We take the Borel mass of 
1 GeV and the self-energies from Ref. [19j . 



In looking for additional sources of significant contri- 
butions at higher densities, we note that as the chemi- 
cal potential increases, the Landau and the unitary cuts 
extend over wider intervals in E, so that contributions 
from intermediate states with higher thresholds need be 
included. Also if the many operators of higher dimensions 
contribute with the same sign, they together may bring in 
a significant contribution. 

It is appropriate here to compare our work with that 
of Ref. [Ij, where similar sum rules are discussed. These 
authors work for p 7^ 0, getting three independent am- 
plitudes and hence three sum rules to determine the two 
self-energies and the nucleon residue. They, however, find 
the latter to be poorly determined by the sum rules. As 
we point out in footnote 1, there is only one combination 
of these amplitudes in the neighbourhood of p = 0, which 
has a large nucleon pole term. We thus write a single sum 
rule for just this amplitude, supplying the (unknown) self- 
energies from independent theoretical determinations. It 
is this focussing on the single large amplitude that makes 
our result for the residue quantitative. 

In our present calculation we face a problem in that 
the coupling parameter A in vacuum itself is not well de- 
termined. This, however, is not the case with the couplings 
of vector and axial-vector currents with the relevant par- 
ticles, though the problem of determining the self-energies 
of these particles in nuclear medium still remains. Once 
the latter problem is treated properly, calculation of such 
parameters may well be quantitative. 

The present work is restricted to nuclear matter at 
zero temperature, so that the pions appearing in the loop 
graphs of Figs. 1 and 2 are virtual. It may, however, be eas- 
ily extended to finite temperature, when there will appear 
real pions, exciting the nuclcons into nucleonic resonances. 



such as Z\(1237) with high degeneracy factors. Then con- 
tributions from intermediate states like 7rZ\ need also be 
included in the spectral side of the sum rule. 

As we work with nuclear matter, we ignore electromag- 
netic interaction. However, in the realistic case of nuclear 
density created, for example, in heavy ion collisions, it 
may not be negligible |36| . 

Finally we put together some similar, known results 
for current-particle couplings and the quark condensate. 
Consider first the pionic medium at low temperature. The 
coupling parameter i^^r of the axial-vector current with 
pion in vacuum changes to F'^ 



12F^) 



(45) 



The coupling of the baryonic current with nucleon, that 
we are considering here, also changes from the vacuum 
value A to [5], 



A^ = A 1 - 



32 F2 I ' 
For the quark condensate, we have 

m 



imV ={0\qq\0) 1- 



(46) 



(47) 



where we keep only the leading term, though it has been 
calculated up to 0{T^) [H]. 

Considering nuclear matter (at zero temperature), the 
Lorentz invariance already breaks at leading order for the 
axial-vector current coupling to pion. 



kf,F„ koF^Sf^o + hF^Sfii 



(48) 



giving rise to two decay parameters. They change with 
nuclear density as [37ll38] 



(49) 
(50) 



F^=.J^^n -(0.26 ±0.04) — 



no 



F^ = F^ll- (1.23 ±0.07) — 

m 



Observe that the two parameters have quite different den- 
sity dependence. But one can argue [37j that it is the tem- 
poral component that reflects the spontaneous breaking of 
chiral symmetry. To these pion decay parameters, we add 
the result of the present work. 



X* = X{l- (0. 20 ±0.04) — 

no 



(51) 



obtained as a linear fit to the curves of Fig. 5 up to normal 
nuclear density. Also the quark condensate in this medium 
is given already by Eq.(37). 

All these results are low-density expansions (in pionic 
and nuclear media) for the current-particle couplings and 
the quark condensate. What is worth noting is that at such 
densities all the couplings along with the quark condensate 
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definitely decrease with rise of density. Tiie significance 
of this decrease is clear. In vacuum, F^, (A) measures the 
overlap of the pion (nucleon) state with those obtained by 
applying the axial (nucleonic) charge on the vacuum. In 
medium their decrease indicates that the external sources, 
A^{x) and ri{x) become gradually weak in exciting pion 
and nucleon in the respective media. The decrease of the 
quark condensate implies a trend towards restoration of 
chiral symmetry. 

It is interesting to recall here that although QCD is 
believed to be the theory of strong interaction amongst 
quarks and gluons, it is asymptotically free at short dis- 
tances. This feature leads to the common belief that with 
the rise of temperature and/or density, hadronic matter 
will undergo one or more phase transitions, restoring the 
spontaneously broken chiral symmetry and liberating colour 
to form quark-gluon plasma. The corresponding order pa- 
rameters are the quark condensate and the Polyakov loop, 
which must be evaluated on the lattice to investigate such 
transitions. Any density expansion of the parameters, even 
if carried to arbitrarily high order, would not be valid close 
to phase transition [55] . 

It is, however, the case that often the above first or- 
der formulae give critical values in qualitative agreement 
with lattice and exact model calculations. The best exam- 
ple is {qq}'^ , for which the value given by Eq. (47) is of 
the same order as the one from the lattice [40]. Thus at 
finite temperature, the same value approximately for the 
coefficients of in Fj, A"^ and (qq)'^ tends to support 
the expectation that they all go to zero at the same (crit- 
ical) temperature [6]. In the same way, at finite nucleon 
chemical potential, we find that the coefficients of n are 
again approximately the same for all the three quantities, 
allowing us to expect that they all disappear together at 
the same critical density, which is several times away from 
the normal nuclear density. Our speculation has an added 
importance in that quantitative calculation on the lattice 
proves difficult at finite chemical potential [H] . 
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